RARE MUTATIONS IN EVOLUTIONARY DYNAMICS 



ANNA LISA AMADORI 1 , ANTONELLA CALZOLARI 2 , 
ROBERTO NATALINI 3 , BARBARA TORTI 2 



Abstract 

In this paper we study the effect of rare mutations, driven by a marked point 
process, on the evolutionary behavior of a population. We derive a Kolmogorov 
equation describing the expected values of the different frequencies and prove some 
rigorous analytical results about their behavior. Finally, in a simple case of two 
different quasispecies, we are able to prove that the rarity of mutations increases the 
survival opportunity of the low fitness species. 
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1. Introduction 

Evolutionary dynamics describes biological systems in terms of three general prin- 
ciples: replication, selection and mutation. Each biological type - a genome or a 
phenotype as well as a species - is described by its reproduction rate, or fitness. In 
force of selection, a population evolves and changes its fitness landscape. Genetic 
changes can help in reaching some local optimum, or open a path to a new fitness 
peak, but sometimes they may dr i ft pop ulation away from a peak, especially if the 
mutation rate is high. See iNowak (|2006l ) for an extensive account of the state of the 



art concerning evolutionary dynamics. 

A mixed population of constant size constituted by a fixed number of different 
types is characterized by the vector collecting the relative abundance of each type: 
x = (xq, x\, . . . , x D ). By virtue of evolution, this vector draws a path in the simplex 

S = < x = (xq,xi ...i D )£l D+1 : x k = 1, x k > as k = 0, 1 . . . D \ . 

I k=0 J 

The mechanism of replication/selection is well described by an ordinary differential 
equation, where the relative fitness measures the balance between death and birth of 
individuals. Denoting by ff. the absolute fitness of any fe-type, by / = (/o, x±,..., f D ) 
the fitness vector, and by / = x ■ f the mean fitness of the population, this equation 
reads 

(1.1) ^- = (fk-f)xk, as k = 0,1,... D. 

Several shapes have been proposed for the absolute fitness. When one is modelling 
phenoty pes, the choice of a constant fitnes s seems fair, but, starting with the seminal 
work by Mavnard Smith and Price] (1973), an impor tant amount of research deal s 



with ideas arising from mathematical game theory, see Hofbauer and Sigmund ( 19981 ) 
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and references therein. In that framework, the individual fitness is taken as a linear 
function of the population x, i.e. f{x) = Ax, where A is the payoff matrix that 
rules the interplay between different strateg ists. In this case, equatio n (jl.ip is the 
celebrated replicator equation introduced by Taylor and Jonker ( 19781 ). 



Concerning mutations, it h as to be mentioned the quasispecies equation introduced 
by Eigen and Schuster ( 19791 ). where constant fitness were considered. This choice 



modifies equation (jl.ip into 

(1.2) ^ =^2fi<lik%i- fxk, as k = 0,1, ...D. 

Here the coefficient express the proportion of offspring of A:-type from a progenitor 
i, which shows up at any procreation. It is clear that q^. = 5ik gives back the equation 
(|1.1|) . When the fitness vector is given by the relation f(x) = Ax, as suggested by 
evolutionary game theory, then equation (jl.2p is the well-kno wn replicator-mutator 



equat ion, also known as selection-mutation equation, studied in lStadler and Schuster 



fll99j). 

As a matter of fact, mutations i ntroduce a rand o m ing redient into evolution, that 
is not enough emphasized in (|1.2|) . Traulsen et al.l ( 20061 ) pointed out that equation 



(jl.2p can be recovered by assuming that the population f ollows a generalized Moran 
proc ess and taking t he lim it for large population size. Champagnat et al. ( 20081 ) 



and iJourdain et al .1 (|2012l ) showed that various macroscopic diffusion models can 
be derived by the same individu al stochastic process, by pe rforming different types 
of rescaling. We also mention iDieckmann and Law where a macroscopic 



dynamics is deduced by an individual based stochastic description of the mutation 
process. 

Here, we prefer to take a more macroscopic viewpoint, which however takes strongly 
into account the different regime of mutation processes. We start by modelling the 
stochastic dynamics at the level of the frequency vector x. In order to capture the 
"rarity" of mutations, we assume that they are driven by some point processes. As a 
result, the random path of the frequency vector x in the simplex is not continuous. 
This happens because mutations arise at a different time-scale with respect to repli- 
cation and selection, and this assumption constitutes the main novelty of the present 
paper. 

The stochastic dynamics for frequencies is introduced and studied in Section [2j In 
Section [3l a Kolmogorov equation describing the expected values is rigorously derived 
and studied in its analytical aspects. Because of the point process, such equation is of 
integro-differential type. Global existence of mild solutions to this equation is proved 
as well as some useful basic estimates. Next, in Section 21 we deal with the particular 
case of two different quasispecies. We prove the existence of a stable equilibrium, 
that in general is greater or equal to the one of the standard quasispecies equation, 
to point out the fact that the rarity of mutations increases the survival opportunity 
of the low fitness species. In this sense, the single relevant contribution of our paper 
is to show that the equilibrium position depends not only on the global amounts of 
mutations, but also on the time intensity of the process driving mutations. As the 
time intensity goes to infinity, the standard quasispecies equilibrium is recovered. 
But when mutation are concentrated in few, very rare events, a different equilibrium 
arises. We prove this occurrence in a rigorous way in Section HI using some refined 
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uniform estimates of the derivatives of the solution, in the simple case of two different 
quasispecies, in a given range of the parameters of the model. 

2. The stochastic model 

We introduce here a stochastic differential equation that describe the population 
distribution, under rare sudden mutations. Let us first revis e the selection-mutation 
equat ion (jl.2p . To underline the effect of mutation, we follow Hofbauer and Sigmundl 



( 19981 ) and introduce the "effective mutation matrix" M = Q — I. If no mutation 



occurs, M is the null matrix; in general M = (raft)i,fc=o,i,...,D with rriik = € [0, 1] 
if i ^ k and mkk = Qkk — 1 = — Yl m ki € [—1,0]. Equation (jl.2p can be written as 

i^k 

d D 

(2-1) —X k = {fk{x) ~ f(x)) X k + ^ fi( X ) m ik x u 

for t > as k = 0, 1, ... D. Here and henceforth we have assumed that all fk(x) 
are nonnegative. We remark that the term (fk(x) — f(x)) x k stands for the homoge- 
neous reproduction steered by the replicator equation, the term ^ fi(x)m i f < .Xi > 

i^k 

describes the increasing of the frequency of fc-type yielded by birth of /c-individuals by 
mutation, and the term fk(x) mkk Xk < measures the decreasing of the frequency of 
k individuals caused by the birth of mutated descendants by fc-type progenitors. The 
underlying assumption is that mutations happen at the same time-scale as homoge- 
neous reproduction: at any procreation, a fixed proportion of the progenies shows 
up a mutant trait. It seems however more realistic to describe mutations as sudden 
changes in the population distribution: we thus assume here that they are driven 
by some marked point processes. To be more precise, let (T n ,Z n ) n be a sequence 
of random times and random marks on a filtered probability space (CI, J-, (J~t), P)- 
The marks are chosen in the mark space Z = {(i,k) G {0, ...,d} 2 , i ^ k}, i.e. it 
is assumed that each mutation has progenitor of a fixed type and descendants of a 
different unique type. Set 

N t = T, I (T n <t), N? = Y.\z n =^))\T n <t), i + h t>0, 

n n 

which define the processes that count respectively the total number of mutations and 
the number of mutations with ancestor i and descendants k, so that Nf = X^fc ^t k ■ 
The intensity of Nl k should depend on the "genetic distance" from type i to type k. 
Besides, it is affected also by selection: the larger the fitness fi(x t -), the more type 
i reproduces, the more often its offspring shall suffer a mutation. Hence we take as 
(J r t)-intensity kernel 

X t (dz) = ^2 Xikfi{x t -)\i,k)(dz), 

i^k 

where Ajfc are positive constants. It follows that, for each i k, N %k is a point process 
with (J^)-intensity equal to Xikfi( x t~)- We remark that when the fitness vector / is 
constant then, for each i ^ k, N is a classical Poisson process of parameter Xikf 1 
(see also Remark 12. ip . 

The proportion of the offspring of ^-individuals showing a fc-type by effect of mu- 
tation is taken constant and denoted by 7^ G (0, 1]. A stochastic dynamics for the 
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relative frequencies arises 

(2.2) dx k = (f k (x t ) - f(x t )) x k dt + J2 7*fc 4dNl k - ]T lki x k dN ki , 

i^k ij^k 

as k = 0, 1 . . . D. We next address to well-posedness of the S.D.E. (|2,2p . In view 
of writing its infinitesimal generator (and, later on, its Kolmogorov equation), we 
introduce the vector valued function a{x) = (ao(x), ai(x) . . .a D (x)), with 

(2.3) a k {x) = - [f k {x) - f{x)) x k , fc = 0,l...D, 
and the first order discrete non-local functional 

(2.4) J(p(x) = ^2 ^ijfi( x ) [4> (x + jijXi(ej - e,)) - cp(x)] . 

Here ej stands for the unit vector pointing in the direction of the j axis. Taking 
into account that for (i, k) ^ (j, I) the processes N lk and N^ 1 have no common jump 
times, Ito's formula for semimartingales gives 

(2.5) d(j){xt) = {-a{x t ) ■ (j){xt) + J<p{x t - )) dt + dM t 

for each function 4> € C^(M D+1 ). Here M is the (J r t)-martingale defined by 

dM t = [<f>(xt- + Hkx\- (e* - e t )) - <t>(x t - )] dM l t k , 

i^k 

with M'i k = N l t k - \ ik j fi(x s )ds. The S.D.E. equation ([22}, endowed with any 

Jo 

random initial position .Fo- measura ble, is well-posed in the weak sense stated by 
next lemma. Moreover it is consistent with frequency modeling, i.e. when the starting 
position is in S then the solution stays in S for all t > 0, a.s.. Actually, the next 
Lemma holds. 

Lemma 2.1. Let f : R D+1 — > R D+1 be a Lipschitz continuous function with compact 
support containing the simplex S and such that \f(x)\ ^ if x G S. Let B be the 
operator defined by 



(x) = -a(x) ■ V<£(x) + ls(x)J(f>(x). 

Then for every probability measure ttq assigned on (R D+1 , £>(K D+1 )), the martingale 
problem for (B,iro) is well-posed, that is there exists a filtered probability space and a 
(d +1)- dimensional Markov process (Xt, t > 0) with initial distribution ttq on it such 
that 

<t>{X t ) - <f>(X ) - [ B^Xs) ds 



is a martingale. This process is unique in law, that is if (Xt, t > 0) is a different 
solution of the martingale problem for (B,no), then X and X have the same finite 
dimensional distributions. 

Moreover, if the support of ttq is contained in S then for all t > 0, Xt € 5 a.s.. 
Proof. For all <f> G C fe 1 (]R D+1 ), I s {x)J(f)(x) can be written as 

X(x) / (f>(x + Is(x)y) - 4>{x) m x (dy) 
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where 

A(x) = } ]\jkfi(x) 
and m x (dy) is the probability measure on (R D+1 , Z3(IR D+1 )) defined by 



(2.6) m x (dy) 



t(x) ^Mek-e^dy), if x G S 
any probability measure, if x g" S 1 . 



Note that by the hypotheses A(x) ^ when x 6 S. 

Moreover it is easy to see by using Skorohod construction for random variables 
that ls{x)J(t){x) can also be expressed in the form 

[4>(x + K(xX))-Hx)} KO 

where (H, T) is a measurable space, (x, £) — > K(x, C) is a measurable bounded func- 
tion on M D+1 x S with values in M D+1 and v(dQ) is a <r-finite measure on (H, T). 

More precisely the previous equality holds for (S, T) = ((0, 1) D x R+, B((0, 1) D ) ® 
£>(K + )) , with general element denoted by Q = (no, «i, ■ ■ ■ u D -i,8), 

v{dQ) = duodui . . . du D _i d6, 

and the function K constructed as follows (see Calzolari and Nappol ( 19961 )). Fixed 
x G S, let Y = (Yq,...,Y d ) be a random vector with law m x (dy) defined by ()2.6p . 
Let yo Fo(yo) be the distribution function of Yq and moreover, for n = 1, . . . , D, 
let F n (y n | j/0) ■ ■ ■ ) 2/n-i) be the distribution function of Y n given Yq = yo, . . . , y n _i = 
y n _i, so that 

m x (dy) = F (dy )F 1 (dy 1 \ y ) . . . F D (dy D \ y , . . . ,y D -i). 

Finally denote by i 7 " 1 the generalized inverse of F n , for n = 0,1,...,D. Then for 
x G S let K{x, ■) : S — > K D+1 be defined by 

K (x,() = F Q - 1 (u ), 

K l (x,C) = F 1 - 1 (u 1 \K (x,0), 

K n (x,() = F- l (u n | Kv{x,C,),Ki{x,Q,...,K n _i(x,Q), n= 1,...,D-1, 

and 

K u (x,0= I(o,A(,.)) W^d -1 (0/\(x) | K (x, C), #i (x, C), • • • , K u _! (x, C)) ■ 

For x ^ S 1 let i^(x, •) be identically zero. 

It is to remark that the above construction implies that, for x G S, we have 

(2.7) i/(C€S,|if(x,O|^0)=A(x). 

In fact the distribution functions used in the construction have no jumps at zero so 
that none of the general inverses yields zero on a set of positive measure. Moreover 
for all x and \K(x, ()\ < 2. Then existence of a solution of the martingale problem 
for (B, ttq) follows by existence of a solution of the SDE 

(2.8) X t =X - [ a(X s )ds+ [ [ K(X s -,()j\f(ds x d(), t>0 

JO JO Jr. 
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where Af(dt x dC,) is a Poisson random measure on (R D+1 , £>(M D+1 )) with mean mea- 
sure dt x y{dC) and Xq is a random variable on the same probability space with dis- 



tribut ion 7Tq. A strong non-explosive Markov solution of (|2.8p exists bv lAthreva et all 



( 19881 ). where more general SDE with jumps are treated. Indeed, under our regularity 



assumption on /, not only the drift coefficient verifies sub- linear growth and Lipschitz 
condition but also the intensity of the point process which counts the total number 
of jumps in (|2.8p is bounded. In fact (|2.7p joint with the regularity of / gives 

supz/(C G S, \K(x, C)| 7^ 0) = sup X(x) < +oo. 

xeS x&S 



Again we refer to Athreva et al. (1 1988!) for deriving uniqueness in law of ( Xt , t > 0) 



Finally, following Athreva et al. ( 19881 ). the construction of the strong solution 



of equation (|2.8p uses sequentially the deterministic and the stochastic part of the 
dynamic. So if the support of ttq is contained in S then every trajectory of the solution 
verifies Xt G S, for all t > 0, when either the deterministic dynamic or the stochastic 
dynamic do not allow to the trajectories starting in S to leave S. As it is well- 
known this is true for the deterministic dynamic. As far as the stochastic dynamic 
is concerned, it is sufficient to note that when x £ S then x + ")ikXi{ek — e%) € S. □ 

Proposition 2.1. The frequencies process is well-defined in law for all t > 0. 

Proof. It follows immediately by previous lemma recalling (|2.5p . i.e. that the fre- 
quencies process starting at x £ S solves for all t > the martingale problem for 
(B,5 X ). □ 

Remark 2.1. Let us remark that when the fitness vector f assumes constant value 
on S, then the point processes N^ k are classical Poisson processes of parameter Xikf 1 
and \2. 2\) can be written in form \2. 8\) with 

5 = Z = {(i, k) 6 {0, ... , d} 2 , i ^ k}, K(x, (i, k)) = I s (x)j ik Xi(e k - e;) 

^([0,t]x(i,fc))=JVf, v({(i,k)})=\ ik f. 
So Athreva et al\ tl98& ) directly applies to our model, and the frequencies process is 



the strong solution of equation 112. 2\) so that 

D 

dt 



dx k t 



{fk(xt) ~ f{xtj) Xt + ^2 fi( X t) X ik7ikX l t 
i=0 

ki 



+ £ llk x\- dMl k - lkix k t - dM\ 

By taking mik = \iklik as i ^ k, nikk = — ^2 ^kilki, the stochastic dynamics ()2.2 

i^k 

is nothing but the quasispecies equation (|2.ip . perturbed by a martingale term. 



3. The Kolmogorov equation 
We next address to the expected value of the frequencies process 

u k (x,t) = £ (x k | x = , 

as k = 0, 1, . . . D, and deduce rigorously that it satisfies its Kolmogorov equation. 



RARE MUTATIONS IN EVOLUTIONARY DYNAMICS 7 

Proposition 3.1. For each k = 0, 1 ... D we have 

d t Uk(x, t) + a(x) ■ Vufc(x, i) = Ju k (x, t), x G S, t > 

u fc (x,0) = x G S 1 , £ = 0, 

We recall that o and J7" have been introduced in (|2.3p and ()2.4p . respectively. 
Proposition 13.11 follows readily by next result about the semigroup on B(R D+1 ) 
defined by 

T t <f>(x) = £(<j>(X?)). 

Here, Xf denotes the solution at time t of (|2.2|) . with deterministic starting position 
x G ]R D+1 . 

Lemma 3.1. Let f and B be as in Lemma \2.1\ Then B is the infinitesimal generator 
of (Tt, t > 0) with domain T>(B) containing C|^(1R £,+1 ). Moreover for each cf> £ 
C\{EP +l ) the scalar function (x,t) — > u ( p(x,t) = T t (j){x) satisfies the Kolmogorov 
equation 

dtu^x, t) + a(x)Vu cj) (x, t) = Ju^x, t), 
u^(x,0) = 4>(x). 

Proof. Any function (j) G C\ belongs to the domain of B since 

T t cj)(x) - 4>(x) 



lim sup 



B(j)(x) 



0. 



In fact Ito's formula joint with Fubini's theorem gives 



(3.1) T t 4>(x) = 4>(x) + T s B4>(x)ds 

Jo 

so that the above limit coincides with 

(3.2) lim t- 1 sup / £ {B(f){X x s ) - B<f>(x)) ds 

*^ 0+ x£R D+1 Jo 



0. 



Then equality (|3.2p follows by considering that for all x G IR D+1 the regularity of / 
and cj) implies that 

S(\B<l,(Xf)-B<l,(x)\)<C(tf>)8, 

with C((j)) a positive consta nt depending on t/>. Finally the thesis follows from (|3.ip 
by recalling that (see, e.g. Lamperti ( 19771 )) if G D(B) then Tt<p G T>{B) and 
BT t (j) = T t B(j), for each t > 0. □ 

3.1. Dimensional reduction. We remark that the number of variable can be re- 

D 

duced by setting xq = 1 — Yl x k- The new variable, that we still denote by x = 

k=l 

(xi . . . x D ), lives in the closed set 

D 

£ = {x = (xi . . . x D ) G R D : x fc > as k = 1 . . . D, ^ x k < 1}. 

fc=i 

In all the following, we continue to write a and / for the respective functions de- 
pending on the D variables x = (xi . . . x D ) G X. As or the non-local term J", it 
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becomes 



J4>(x) = ^2 Xijfi(x) {(f) (x + 7ijXi(ej - ei)) - 4>(x)} 
ij'=l...D 

+ XI X iofi(x) [<j) (x - Tio^ej) - <j){x)} 



i=l... D 



j=i...d 



V i=l...D / 



4>(x) 



for any continuous scalar function <p G C(E;K). 

Next subsections are devoted to the analytical study of the decoupled system 

dtUk(x, t) + a(x) ■ Vu k {x, t) = J'ukix, t), x G X, t > 

u fc (x,0) = Xfe, x € S, i = 0, 



(3.3) 
as k = 1 . . . d 



3.2. Mild solutions. We establish global well-posedness of the Cauchy problem 
()3.3p . As a preliminary, we notice that the vector field a(x) does not point outward 
at the boundary of S. 

Remark 3.1. The projection of the vector field a orthogonal to the sides of the 
boundary o/S is always zero. Actually at x k = we have a(x) ek = a k (x) = 0, while 

D 

at ^2 Xk = 1 we have 
fc=i 

(D \ D D D 

X e k j = X a k (x) = -^x/i fk(x) + f{x) ^x k = -xf{x) + f(x) = 0, 
k=l J k=l k=l k=l 

D 

because xo = 1 — ^ x k = 0. 

fe=i 

As a consequence, the flux of the Cauchy problem for the autonomous equation 

(3.4) { y = a ^ 

1 2/(0) = x, 

is well defined and maps S into itself. It is worst noticing that it is nothing than 
the solution to the replicator equation (jl.ip . after the dimensional reduction. In the 
following, we shall write Y(x, t) for the solution of (j3.4|) starting at x. For any x, t, 
the function s t— > Y(x,s — t) is the characteristic line of through x, t for problem (|3.3p . 

Definition 3.1. A mild solution is a function u G C ([0, T) x £;]R D ) satisfying the 
integral formula 

(3.5) u k (x,t) =Y k (x,-t) + Ju k (Y(x,s -t),s)ds, 

Jo 

as k = 1 ... D. 

Following the method of characteristics, we define 

v(x, s, t) = u(Y(x, s — t),s) 
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(3.6) 



and notice that for any t, v solves an integro-differential problem 

d s v k (x, s, t) = lv k (x, s,t) x G S, s > 0, 

v k (x,0,t) = Y k (x, —t) x G £, s = 0, 

as k = 1 ... D. Here 

lv k (x,s,t) = ^ ^ijfi( Y ( x , s - t)) [v k (x + 5ij(x,s - t),s,t) - v k (x,s,t)] 

i,j=l...D 

+ ^2 X i0fi( Y ( X ' S ~ *)) [ V k(x + $io(x,S - t),S,t) - V k (x,S,t)] 
i=l...D 

+ X] ^0jfj(Y(x, s - t)) [v k (x + 6 0j (x,s - t),s,t) - v k (x,s,t)] , 
j=i...d 

Sij(x,r) =Y (Y(x,r) + j ij Y i (x,r)(e i - e ; ) , -r) - x, 
5 i0 (x, r) =Y (Y(x, r) - j i0 Yi(x, r)e ; ) , -r) - x, 

S j0 (x,r) =Y(Y(x,r) + 7oj(l - ^ r ^ e J ' _r ) ~ x ' 

1=1... D 

The Cauchy problem (j3.6j) is seen in the mild sense, actually 

(3.7) Vk(x,s,t) =Y k (x,-t) + / Iv k (x,r,t)dr. 

Jo 

We next use a fixed point argument to solve (|3.6p . This brings a mild solution also 
to (|3.3p because 

v k (x,t,t) = u k (x,t). 
The following proof is standard, but we report it for the sake of completeness. 

Proposition 3.2. The problem (|3.3p has an unique global solution u G C (S x [0, +oo) 

Proof. Let 5, T > 0, and % the set of continuous scalar functions of (x, s,t) G £ x 
[0, 5] x [0, T], which is a Banach space endowed with the sup- norm. For k = 1 . . . D 
and r > , let be the set of functions of x such that 

Bixp{\\v(x,8,t) -Y k (x,-t)\\ : x G S,0 < s < 5,0 < t < T} < r, 

and T the operator 

Tv(x,s,t) = Y(x, —t) + / Xv(x,r,t) dr. 

Jo 

T maps B r into itself because for v £ B r we have 

HTu - foil < 211/HooslbH < 2||/|| 00 ,S'(r + ||u ||) < r, 

provided that 5 < r/2||/|| tX) (r + ||u ||)- Moreover T is a contraction since 

r 

\\Tv — Tw\\ < 2\\f\\ 00 s\\v — w\\ < 2 / WoqS \\v — w\\ < n — rr w — toll. 



r + \\v \ 

It follows by contraction Theorem that (|3.6p has an unique mild solution u G 
C([0,S] x [0,5] x E) at least for 5 = r/2\\f\\ OQ (r + ||u ||). Moreover, since v G B r , 
we have ||u|| < r + ||i> ||. Then the fixed point argument can be iterated to get, 

at any step n, a solution defined until 5 n = (^ r+ \\ Vo \\ + Yl (fc+i)"J+ || l || J ■ As 

S n — > +co, existence and uniqueness of a global solution follows. □ 
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3.3. Viscosity solution. Using the tools of viscosity theory, one can eventually 
prove the following result 

Proposition 3.3. Let u and v be, respectively, viscosity sub/ supers olutions of (|3.3p 
with u(x, 0) < v(x, 0) for all x € X. Then u(x, t) < v(x, t) for all x € £ and t > 0. 

Proof. The proof consists in assuming that there exists T > so that 

M = max{e _i (tt(x, t) - v(x, t)) : < x < 1, < t < T} > 

and getting a contradiction. Standard arguments exclude that any maximum point 
x, i may have i = either i € (0, T] and x in the interior of X. We show that neither 
x in the boundary of X is allowed. To fix idea, suppose that x is a strict maximum 

D 

point with x\ = 0, Xi > as i = 2, . . . D and ^ Xj = 1. Next, a barrier function 

i=i 

D 

k(x, t) = 1/xi + 1/(1 - ^ x0 + 1/(T - t) 

i=l 

is introduced and M is approximated by M(5), the maximum value of 

e~ l (u(x, t) — v(x, t)) — 8k(x, t). 

It is clear that there exists a maximum point x(5),t(5) with x(S) in the interior of S 
and < t(5) < T. Moreover 

M(5)^M>0, x(5)^x, t(d)^i>0, 6 n(x{5),t(5)) ->■ 

as S — > 0. Now the perturbed functions u$(x, t) = e~'n(x, t) — 5 k(x, t) and ^(x, t) = 
e~ t v(x,t) can be handled by the standard tool of doubling variables and using equa- 
tion (|3.3p . The step of passing to the limit as 5 — > can be performed since 

/ 



\5a(x(6)) -D x K(x(5),t(5))\ < 



ai(x(S)) 



\ 



xi(5) 



+ 



1=1 



Y,(xi - Xi(5)) 



i=i 



5K(x(S),t(5)), 



where the quantities \ai(x(6))/xi(5)\ and \Y1 a i( x ($)) / Yl( x i ~ x i($))\ are bounded 



i=l 



i=l 



because a is Lipschitz continuous and a\{x), ^ a,i(x) are zero by Remark 13.11 □ 



8=1 



In our particular setting, we have decided to emphasize the transport component 
by following the method of characteristics. Nevertheless, the solution defined and 
produced in the previous section coincides with the viscosity solution. 

Proposition 3.4. Mild and viscosity solution of (|3.3p coincide. 

Proof. As well posedness holds in both framework, it suffices to check that the mild 
solution is a solution in viscosity sense. We only prove the subsolution part, because 
the supersolution is identical. Let k = 1 . . . D, and <f> be a smooth scalar function such 
that Uk — 4> has a global maximum at (x,t) G E x (0, +oo), with Uk(x,t) = 4>(x,t). 
First, we notice that Ju^{x, t) < J<fi{x, t). Next, we set ip(x, s, t) = <fi(Y(x } s — t), s): 
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it is clear that also v k — ip has a global maximum point at (x,t,t), and therefore 
Ivk(x,t,t) < Xip(x,t,t). Hence it is easily seen that 

d t (f>(x, t) + a(x) ■ D x <j>(x, t) = d s ip(x, t, t) = d s v k (x, t, t) 

= lv k (x, t, t) = Ju k (x, t) < J(p(x, t). 

□ 

Comparison techniques provides furth er information about the regularity of u 
w.r.t. x (see, for instance, Amadori ( 20031 )). 

Proposition 3.5. The solution to (|3.3|) is Lipschitz continuous w.r.t. x, for every 
fixed t > 0. 

Remark 3.2. It is not hard to extend all results in this section to a more general 
class of problems that can be written as (|3.3p . where J is a first order non local 
operator in the form 

j<p{x) = f f(x,C) [4>(x + 1 (x,c))-4>(x)} du(0, 

provided that 

i) £ is a compact set, whose boundary is globally Lipschitz, and consists in the 
union of smooth surfaces which have exterior normal vector. 

ii) a is a Lipschitz continuous vector field defined on S. At every x in the boundary 
of S, and for every exterior normal vector n, we have a ■ n = 0. 

Hi) v is a finite measure. 

iv) j is a (vector valued) function, Lipschitz continuous w.r.t. x G T, (uniformly 

w.r.t. C) with x + 7(x, £) G £ for all x S S and v-almost any £. 
v) f is a (real valued) function, Lipschitz continuous w.r.t. x € £ (uniformly 
w.r.t. () with f(x, C) > for all x € £ and v-almost any (. 

vi) u is a continuous (vector valued) function. 

4. TWO QUASISPECIES 

With the aim of expounding the behavior of the expected frequencies, we deal with 
the simplest case: two species (i.e. D = 1) and constant fitness (i.e. fi(x) = fi, as 
i = 0, 1). We also introduce the selection rate 

s = fo - h- 

To fix ideas we take S > 0, so that Xq = 1, x\ = is the only asymptotically stable 
rest point for the replicator equation (jl.ip . 

After the dimensional reduction x = x±, xq = 1 — x, the (scalar) Kolmogorov 
equation (|3.3p for 

u(x,t) = £ (x(t) | x(0) = x) 

reads 

{d t u + s(l - x)xd x u = XofoJou + XifiJiu, < x < 1, t > 0, 
u(x,0) = x 0<x<l., 
where 

Jqu{x, t) = u{x + 7o(l — x),t) — u(x, t), J\u{x, t) = u(x — 71X, t) — u{x, t). 

Here and henceforth we have shortened notations by writing 70 = 701, 71 = 710 1 
Ao = A01, Ai = Aio- 



(4.1) 
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Our main concern is to compare the dynamics with rare mutation and the stan- 
dard quasispecies dynamics. In this particular setting (and after the dimensional 
reduction) the quasispecies equation (|2.ip reads 



(4.2) 



x = — S x(l — x) + mo/o(l — x) — mi fix t > 0, 
x(0) = x. 



Here mo stands for the mutation parameter from species to species 1, and viceversa 
for mi. With Remark 12.11 in mind, we take 

m% = Aj7j, as i = 0, 1. 

Let X(x,t) stand for the flux associated to (|4.2p . i.e. the solution to the o.d.e. with 
initial condition x, computed at t. It solves the homogeneous transport equation 



(4.3) 



d t X + (ss(l-i) - m f (l-x) + mi fix) d x X = 0, < x < 1, t > 0, 
X(x,0)=x, < x < 1. 



Let us recollect some basic facts about (|4,ip . In this simple scalar setting, u is a 
classical smooth solution. 

Proposition 4.1. The Cauchy problem (|4.ip admits a classical solution u G C°°([0, 1] x 
[0,oo)). 

We do not report in details the proof of this result because it is completely standard: 
it relies in deriving w.r.t. x iteratively the equation and the initial datum in (|4.ip and 
noticing that the obtained problem inherits the same structure and regularity. We 
rather go into details and obtain some more estimates concerning first and second 
order derivative w.r.t. x. 

Lemma 4.1. For every t > 0, the solution to (|4.ip satisfies 

< d x u(x,t) < e ( s - m o/o-mi/i)< 

for all x € [0, 1]. 

Proof. Let p = d x u: deriving (|4.ip w.r.t. x gives 

{d t p + s(l - x)xd x p+ (s(l - 2x) + m f + mifi)p 
= A /o(l - To)JoP + Ai/i(l - 7i)JiP, 
p{x,0) = 1. 

It is clear that p = and p = e^ s - -m o/o-mi/i)i are ^ respectively, a subsolution and a 
supersolution. So the thesis follows by comparison. □ 

A similar estimate holds also in the general case treated in previous section, and 
can be proved for viscosity solutions. It is worst mentioning that, in particular, u is 
monotone increasing w.r.t. x, for every fixed t > 0. It is also convex, as shown by 
next lemma. 

Lemma 4.2. For every t > 0, the solution to (|4.ip is convex w.r.t. x. Moreover 
there exist two constants C > and ^ 6 R such that 

Q<dl x u{x,t)<ce^ 

for all x € [0, 1] . 
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Proof. Deriving (|4.ip twice w.r.t. x gives that q = d xx u solves 

d t q + S x{l — x)d x q + (4 s(l - x) + a)q = 2 S d x u+ 

Ao/o(l - To) 2 Joq + Ai/i(l - "fi fjiq, 

with a = —2s +(2 — 7o)w,o/o + (2 — 71)^1/1 ■ As c^u > 0, the function d xx u is a 
supersolution to the homogeneous Cauchy problem 

f d t q + Sx(l- x)<9 x g + (4S(1 - x) + a)q = A /o(l - lo?Joq + Ai/i(l - lifJiq, 
1 ?M) = 0, 

and therefore d xx u > 0. On the other hand, as d x u < e^ s - m o/o-mi/i)t^ ^ e f unc ^ n 
e?;Lu is a subsolution to 



(4.5) 



d t q + Sx(l - x)d x q + (4 s(l - x) + a)q = 2 S e^ 8 -"»o/o-mi/i)t + 

Ao/o(l - 7o) 2 Jo? + Ai/i(l - nfJiq, 

j(x,0) = 0. 



Eventually, also the estimate from above of follows by comparison, after having 

2 S 

checked that the function q(x,t) = — e (2s -m o f -m 1 f 1 +e)t - g a SU p erso lution to (j4.5|) . 
for every e > 0, □ 

It follows that the expected value of the density with rare mutations is greater or 
equal than the deterministic one, i.e. rare mutations increase the survival opportuni- 
ties of the low-fitness species. 

Proposition 4.2. Let u and X be, respectively, the solution to (|4.ip and (|4.3|) . Then 
1 > u(x, t) > X(x, t) for all < x < 1 andt>0. 

Proof. It is clear that u(x,t) < 1, because the constant function 1 is a supersolution 
to (|4.ip . Concerning the estimate from below, it follows by comparison after checking 
that it is a supersolution of ()4.3p . Indeed 

d t u + (sx(l - x) - m /o(l -x) + mi fix) d x u = 
Ao/o [u(x + j (l - x),t) - u(x,t) - 70(1 - x)d x u(x,t)] 
+A1/1 [u(s - jix, t) - u(x, t) + 7ix9 x -u(x, t)] = 

Ao/o / d xx u(x + 6*70(1 - x),t)d9 + Ai/i / d xx u(x - Ojix, t)d0 > 
Jo io 

by convexity. □ 

4.1. Large time behavior. It is well known that the quasispecies equation (|4.2p 
has an asymptotic equilibrium at the point x € [0, 1] singled out by the relation 

S x(l — x) — mo/o(l — x) + mifix = 0, 

and that its basin of attraction is given by [0, 1] or [0, 1), depending on the value of the 
parameters (see Remark 14.11 later on). In our notation, this means that the solution 
to (|4.3p satisfies lim X(x,t) = x for every x 6 [0,1] (or for every x 6 [0,1)). 

t— >+oo 

It is interesting to study whether the rare mutation equation (|4.ip has the same 
asymptotic behavior, or rather exhibits a new equilibrium. To begin with, we need 
to establish that the solution to (|4.ip actually admits a limit as t — > +oo. For a 
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restricted range of parameters, the estimate given in Lemma 14.11 suffices to deduce 
the large time behavior of u. Otherwise, some more work is needed. 

Lemma 4.3. If s > r?io/o + mifi, then the solution to (|4.ip satisfies 

d x u(x,t) < -. FTTe- mi/lt 

y ' ' - 1 — x + e -(s -m fo)t 

for all x € [0, 1] and t > 0. 

Proof. It suffices to check that p = 2e ( - s - m ° fo+mif ^ t /(l + (1 - x)e( s " m o/o)<) is a 
supersolution to (|4.4p . To shorten notation, we write z = (1 — x)e ( - s ~ m °f°^ t . Trivially 
0) = 2/(2 — x) > 1. Moreover easy computations give that 

«9 t p + s(l - x)x d x p + (s(l - 2x) + m f + mifi)p 
9 e (s -mo/o+fni/i)i 

-(2s(l-a0 + ((s(l-x) + mo/o)z) 



^ , .a 2 m o/o^, 



(1 + z) 2 
2 e (s -m /o+mi/i)i 

(1 + z 

JlP < o, 

2 e (s -mo/o+"n/i)t 2e^ s _m o/o+nii/i)t 

JbP = — , r-, rr7o^ < 



(l + z)(l + (l- 7o )z) " (1 + z) 2 l- 7o ' 
and the thesis follows immediately. □ 

We are now in the position to draw the large time behavior of u. 
Proposition 4.3. For every choice of the parameters, the function 

u(x) = lim u(x, t) 

t—>+oo 

is well defined for every x G [0, 1]. //, in addition, mi > or mi = and mo > S / fo, 
then the limit u is constant. 

Proof. If mo/o + Tii/i > S, it follows by Lemma 14.11 via standard arguments that u 
converges to a constant as t — > +oo (uniformly w.r.t. x). 

Similar statement follows by Lemma [4.31 if mo/o + m ifi < S and mi > 0, provided 
that x stay in any closed subset of [0, 1). Concerning the behavior at x = 1, we 
deduce by equation (|4.1|) that 

d t u(l, t) = Xifi (u(l - 71, t) - u(l, t)) < 0, 

as u is increasing w.r.t. x. Hence u(l) = lim u(l,t) exists and is finite, actually 

t— >+oo 

r+oo 

tt(l) = 1 - Ai/i / [u{l,t)-u{l-ji,t)}dt. 
Jo 

In particular 

lim (u(l,t) - u(l - 71, i)) = u(l) - u{\ - 71) = 0, 

i— >+oo 

because the function t i— > — u(l — 71, t) has limit as t — > +00 and has finite 

integral on [0, +00). This, in turns, implies that u is constant up to x = 1. 
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For the case uiq = s //o and mi = 0, we know by Proposition 14.21 that 1 > u(x, t) > 
A(x,i) for every x G [0,1] and i > 0. But in this special case the asymptotically 
stable point of (|4.3|) is sc = 1, hence 

1 - x 



1 > u(x,t) > X(x,t) = 1 



l + (l-x)si' 
and u = 1. 

The proof is completed by checking that lim u(x, t) exists even if mo < S //o and 

t— >+oo 

mi = 0. In that case, it follows by Lemma 14.31 that u(x,t) is equicontinuous w.r.t. x 
in any closed subset of [0, 1). Thus standard machinery for evolution equations yields 
that u is equicontinuous w.r.t. both x and t and therefore u(x) is well defined (and 
continuous) for x G [0, 1). On the other hand equation (|4.ip states that dtu(l, t) = 0, 
so that u(l) = 1 = u(l,t) for all t. □ 



Remark 4.1. We mention in passing that Lemmas \4-1\ and \4.3\ imply that u(x, t) —> u 
uniformly w.r.t. x G [0, 1] if S < mo/o + mifi, or, respectively, uniformly w.r.t. x 
in any closed set contained in [0, 1), if S > mo/o + f^xfi- This behavior reflects that 
one of X(x,t) —> x. It is worth noting that, in case mo < S //o and mi = 0, the 
quasispecies equation (j4.3|) gives 



X(x,t) = x + 



TOi/i x = mo/o / S < 1. Therefore the basin of attraction of x is only the interval [0, 1) 
and even the asymptotic limit of X jumps from x to 1 at x = 1 . 

For some choice of parameters, rare mutations give the same equilibrium of con- 
tinuous mutations. This happens, for instance, if mo = 0. In this case the mutated 
descendants have higher fitness than their progenitors, and mutation helps selection 
in fixing the high-fitness specie. 

Proposition 4.4 (Fair mutation). Assume that mo = 0, so that the equilibrium for 
both the quasispecies dynamics (|4,3p and the replicator equation (jl.ip is x = 0. The 
same holds also for (|4.1|) . i.e. u = 0. To be specific, we have 

X e- st > u(x,t) > e -(B+m l/l)t 

1 - x(l - e bt ) 1 _ e -(s+mi/i)) 

S+mi/i 

for all t > 0. 

It has to be remarked that the first and last terms of the inequality are the solution 
to the replicator equation (jl.ip . and quasispecies equation (|4.3p . respectively. 

Proof. As mo = 0, the Kolmogorov equation (|4.ip becomes 

dtu + s(l — x)x d x u = Ai/i [n((l — 7i)x, t) — u(x, t)] < 

because u is increasing w.r.t. x. Then u is a subsolution of the transport equation 
dtu + s(l — x)xd x u = and the first inequality follows. In particular, we have that 
lim u(x, t) = for all < x < 1. As for x = 1, we know that 

t— s>+oo 

u(l, t) = 1 - Ai/i / (tt(l, s) - u(l - 71, s)) ds. 
J o 
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Since the integrand is nonnegative by monotonicity, the function t 1— > u(l, t) is mono- 
tone decreasing and bounded, so it converges. Hence the function t i— >• u(l,t) — u(l — 
71, t) is nonnegative, has finite integral in [0, +00) and has limit as t — > +00. Even- 
tually lim u(l,t) = lim u(l — 71, i) = 0. The proof is now complete because the 

t— >+oo t— >+oo 

second inequality has been established in Proposition 14.21 □ 

The large time behavior of rare mutations reflects the one of continuous mutation 
also in the opposite situation, i.e. when mutation towards the low-fitness specie is so 
relevant to overwhelm selection. 

Proposition 4.5 (Unfair mutation, strong case). Assume that mi = and mo > 
S //o; so that the equilibrium of the quasispecies dynamics (|4. 3[) is x = 1. The same 
holds also for (|4.ip , namely = 1. Moreover 

1 — x 

l>u(x,t)>l — - — «/m >s// , 

1 + °LL x > ( e mo/ -s t _ 1) 

or 

l>"(x,t)>l- 1 + ^; g)t <fmo = s// . 

Notice that the quantity in the right-hand side of both inequalities is the solution 
of the respective quasispecies equation. 

Proof The thesis follows by Proposition 14.21 because in this particular setting the 
equilibrium condition for the standard quasispecies equation reads (sx-mo/o)(l-x), 
and the only root contained in the segment line [0, 1] is x = 1. □ 

Something new happens when mutation is unfair (i.e. m\ = 0) but too weak to 
overwhelm selection (i.e. < mo < s //o)- In this case, the behavior at large time 
depends on the time intensity of the point process governing mutations, and it does 
not follow the relative quasispecies equation anymore. As expected, the quasispecies 
equation is recovered as the time intensity goes to infinity. This topic is illustrated 
in next subsection. 

4.2. Unfair, but weak, mutation. We go into more details and inspect the case 
mi = 0, S > rno/o- As only the coefficients /o, mo, 70 and Ao have effects, we shall 
omit to write the index "0". The quasispecies equation (|4.3p reads 

(4.6) d t u + {sx-mf)(l-x)d x u = 0, 

and has a stable rest point at x = mf / s. Its solution can be explicitly written as 



^ mf 

X(x,t) = ^ + 



mf ( 



S l + Mlz4( e (s-m/)t_D' 
s — mf v 1 

and only depends by the parameters S and mf. In the rare mutation setting, there is 
an entire curve of parameters (7, A) € (0, 1) x (m, 00) that give back the same m and 
S: this curve can be seen as the graph A = m/7. As 7 goes to 0, the time intensity A 
increases, and the paths of the point process driving mutations becomes continuous. 
On the contrary, at 7 = 1 the time intensity gets its minimum A = m, and mutations 
are concentrated in rare events that happen simultaneously to all individuals. The 
respective Kolmogorov equation is 

(4.7) d t u + s x(l - x)d x u = mf [u(l, t) - u(x, t)] . 
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It is easy to check that the only solution with u(x, 0) = x is 

(1 _ x \ e (s-mf)t 

ZM) = 1 "l + (l-x)( e ^-l)- 

In order to study the dependence of the expected density by the parameter 7 (equiv- 
alently, by the time intensity A = m/j), we denote by u 7 the respective solution of 
(|4.ip . namely 



(4. 



TTlf 

dfUj + Sx(l — x)d x u 1 = c7om 7 , < x < 1, t > 

■u 7 (x,0) = x, < x < 1. 



The graph of (0, 1] x [0, 1] x [0, +00) B (7, x, t) h-> ii 7 (x, t) is a continuous hypersurface 
that spans the region between the graph of X(x,t) and the one of Z(x,t). 

Proposition 4.6. For every (x,t), the function (0, 1] 3 7 h-> -u 7 (x,t) is nondecreasing 
and continuous, with ui(x,t) = Z(x,t) and limu 7 (x,i) = X(x,t). Moreover both 

continuity and convergence are uniform w.r.t. (x,t) in each compact set [0, 1] x [0, 7]. 

Proof. To begin with, we check that the functions u 7 are continuous and ordered 
w.r.t. 7. A (formal) derivation of equation (|4.8|) yields that w(x,t,'f) = d^ 1 u 1 {x,t) 
solves 

m f 

d t w + Sx(l -x)d x w = —Jow + h 0<x<l,t>0,0<7<l, 
(4.y) < 7 



w(x, 0, 7) = < x < 1, t = 0, < 7 < 1, 



where 



/i(x, t, 7) =^2" [ u 7( x > *) ~ u -y( x + 7(1 ~ + 7(1 ~ x)d x u~/(x + 7(1 - x), i)] 



^(1-x) 2 f ed^ x u 7 (x + e^(i - x ),t)do. 

* Jo 



'0 

By Lemma 14.21 < /i < Ce' 4 * (with, possibly, a different constant c). Hence 
comparison principle gives < w < Ce^*. This yields, in turn, that the function 
(0,1] 971—)- u*y(x,t) is nondecreasing and Lipschitz continuous, and furnishes an 
estimate of the Lipschitz constant, which is equibounded for all (x, t) £ [0, 1] x [0, 7], 
as 7 > 0. In particular, u 7 gets near Z as 7 — > 1, with uniform convergence for 
(x,t) G [0,1] x [0,71. 

We next check that it 7 approaches X(x,t), as 7 — >• 0, by the viscosity solution ap- 
proach. Let us begin by defining 

u + (x,t)= limsup Uy(y,s) and u~(x,t)= liminf u 7 (y,s). 

(y,s,7)-Kx,t,o) " (y,s,7)-Ka,*,o) 

By construction, u + and u~ are respectively upper and lower semicontinuous, more- 
over u + (x,t) > u~(x,t), and certainly ii + (x,0) = x = u~(x,0). It is trivial to 
check that u + and u~ are (possibly discontinuous) viscosity sub and supersolu- 
tion to the transport equation (|4.6p . Therefore by comparison u + < u~. Thus 
u + (x,t) = u~(x,t) is continuous and equal to lim u 7 (x,t). Next, uniqueness for the 

7->0 

transport equation yields that limt( 7 (x,i) = X(x,t) pointwise. Eventually, Dini's 

7->0 

monotone convergence Theorem implies uniform convergence on any compact set 
[0,1] x [0,7]. □ 



18 A.L. AMADORI, A. CALZOLARI, R. NATALINI. B. TORTI 

The family Uy(x,t) spans the segment between X(x,t) and Z(x,t). As time in- 
creases, the quasispecies solution X(x,t) converges to the equilibrium point x = 
mf I S < 1, while Z(x,t) — > 1. Similarly we expect that the asymptotic equilibrium 
of u 7 spans the segment between x and 1, as 7 goes from to 1. To this aim we 
investigate the large time behavior of the functions ii 7 . Trivially lim ii 7 (l,i) = 1 

t— >+oo 

for any 7. Besides lim u 7 (x,i) does not depends by x 6 [0, 1), for all values of 7 

t— !-+oo 

except at most one. 

Proposition 4.7. Take S > mf and let 7* € (0, 1) be the only solution to 

S 7 + mf log(l — 7) = 0. 
If j S (0, 1) \ {7*}, then there exist a number u 7 so that lim n 7 (x, t) = u~ f for every 

t— >+oo 

x e [0,1). 

Proof. We establish that for every 7 G (0,1), 7 / 7*, there exist a £ (0, 2) and (3 > 
such that the solution to (|4.ip satisfies 

(4.10) d x u(x, t) < 2e- /3t / ((1 - x) a + e" (s " m /+^ 

The thesis follows by (|4.10p by the same arguments of Proposition 14.31 In view of 
proving ()4.10p . we follow the line of Lemma 14.31 and check that, for any a G (0,2), 
there exists /3(a) £ M (possibly negative) such that 

p = 2e( s " m /)'/(l + (1 - x ) a e {s -nf+W) 

is a supersolution to (|4.4|) . Set z = (1 — x) a e^ +s ~ m f) t 7 we have by computations 

dtp + s(l — x)xd x p + (s(l — 2x) + mf)p 
2e (s -mf)t 



2ze (s-mf)t 



2 s(l - x) - s(2 - a)x + (s +m/ - 0) z) 



Z , ^ 2 (s(o - 1) + mf - p) , 



2z e( s ~ m -^)* 
(1 + z) 

(1-(1- 7 )")(1+^) Q . . . . 

tor k(z = ; : . hmce k is monotone increasing we get 

v; 1 + (1 -7)«z 

2„ p (s -m/)i 

(1 - 7) 1_a - 1 

Hence p is a supersolution provided that j3 < /3(a) = S(a — 1) — m/ . 

7 

We conclude the proof by showing that [0, 2] 3 a 1— > /3(a) has a positive maximum. 
Indeed, /3(a) is strictly convex with /3(1) = 0, therefore its maximum is positive 
unless it is reached at a = 1. But a = 1 is not a critical point for 7 / 7*, because 
/3'(1) = s + log(l -7)m//7 / 0. □ 

We already know that mf j S < u 7 < 1; actually we can prove more, namely that 
n 7 — > mf j S as 7 — > 0, and ?2 7 — > 1 as 7 — )• 1. 
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Proposition 4.8. We have lim u y = 1, and lim u 7 = Indeed, for every e > 

■y-¥l 7-5>0 b 

and L G (0, 1), there are T > and To, I\ € (0, 1) so that 

l-s< u 7 (x,t) < 1 for all 7 G [r^l], x € [0,1], t > T, 

x < u 7 (x,t) < x + e for all 7 G (0,r ], x G [0,L], t > T. 

Proof. We first deal with 7 near 1. The function (x,i,7) 1— >• u 7 (x,i) is monotone 
increasing both w.r.t. x and 7. Therefore u 7 (x,i) > ur(0, i) for all 2 G [0,1] and 
7 G [r, 1]. Besides also 1 1-> ur(0, t) is monotone increasing w.r.t. t because by (|3.5|) 

u r (0,t) = -f / (n r (r,s) -n r (0,s))ds 
1 Jo 

with nr(r,s) > nr(0,s) for any s. Hence u^(x,t) > nr(0,T) for all x G [0,1], t>T 
and 7 G [r, 1] and the statement is proved by exhibiting T and T so that ur(0, T) > 
1 — £. But, since ui — > 1 as t — > +00, there exists T such that u\(0,T) > 1 — e/2. 
Next Proposition 14.61 ensures that there is T so that ttr(0, T) > ui(0, T) — e/2 > 1 — e 
as desired. 

Concerning the behavior for small 7, we may assume without loss of generality that 
x < 1 — e. We next perturb the mutation coefficient by means of m e = m(l + es /2/), 
and denote by X £ the solution of the corresponding quasispecies equation (|4.6p . It is 
easily seen that 

X £ (x,t) < mf/s+e 

for all x G [0, L] and t >T, provided that we chose T sufficiently large. 

Thus the thesis follows by comparison, if we exhibit T such that X e is a superso- 
lution to (USD for any 7 G (0,T]. But 

d t X £ + s x(l - x)<9 x X £ - [X £ (x + 7(1 - x), t) - X £ ] 

7 

ems('l + (l-(l + 2|)7>) 
2(l + z) 2 (l + (l- 7 )z) 

where we have used the notation z = fe( s_me ^^ — l) > 0. Taking T = 

S-m £ f v 1 - 

2/ 

1/(1 H ) ends the proof. □ 

£ s 

Eventually, if the point process driving mutation has small intensity (i.e. if 7 is near 
1), the relative population density does not tend to the quasispecies equilibrium as 
t — > +00: the asymptotically stable strategy according to the quasispecies equation 
is not asymptotically stable in expectation, according to rare mutation. 
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